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ABSTRACT 

We present our study on the infrared variability of point sources in the 
Small Magellanic Cloud (SMC). We use the data from the Spitzer Space Tele¬ 
scope Legacy Program “Surveying the Agents of Galaxy Evolution in the 
Tidally Stripped, Low Metallicity Small Magellanic Cloud” (SAGE-SMC) and 
the “Spitzer Survey of the Small Magellanic Cloud” (S 3 MC) survey, over three 
different epochs, separated by several months to three years. Variability in the 
thermal infrared is identified using a combination of Spitzer *s IRAC 3.6, 4.5, 5.8, 
and 8.0 /irn bands, and the MIPS 24 /irn band. An error-weighted flux differ¬ 
ence between each pair of three epochs (“variability index”) is used to assess 
the variability of each source. A visual source inspection is used to validate the 
photometry and image quality. Out of ~2 million sources in the SAGE-SMC 
catalog, 814 meet our variability criteria. We matched the list of variable star 
candidates to the catalogs of SMC sources classified with other methods, avail¬ 
able in the literature. Carbon-rich Asymptotic Giant Branch (AGB) stars make 
up the majority (61%) of our variable sources, with about a third of all of our 
sources being classified as extreme AGB stars. We find a small, but significant 
population of oxygen-rich AGB (8.6%), Red Supergiant (2.8%), and Red Giant 
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Branch (<1%) stars. Other matches to the literature include Cepheid variable 
stars (8.6%), early-type stars (2.8%), young-stellar objects (5.8%), and back¬ 
ground galaxies (1.2%). We found a candidate OH maser star, SSTISAGE1C 
J005212.88-730852.8, which is a variable O-rich AGB star, and would be the first 
OH/IR star in the SMC, if confirmed. We measured the infrared variability of a 
rare RV Tau variable (a post-AGB star) that has recently left the AGB phase. 
Fifty nine variable stars from our list remain unclassified. 

Subject headings: galaxies: Magellanic Clouds - infrared radiation - stars: vari¬ 
ables: other stars: AGB and post-AGB - stars: formation - stars: early-type 


1. Introduction 

Variability is a common feature of many astronomical objects, among which are pulsat¬ 
ing stars, disks around Young-Stellar Objects (YSOs), accreting binary systems, and Active 
Galactic Nuclei (AGNs). The processes behind variability phenomenon are di verse: t h ey 


range fro m sporadic changes in the accretion rate of YSO and AGN disks (see iFlahertv et al. 


20121 and iPeterson et al.l |2000L respectively) to the periodic variat i ons of the opacity and 


2011 


ionization states in pulsating stars’ interiors fllta fe Matsunaeal 12011c iMatsunaga et al. 
Mad ore et al.l 2009). Monitoring the brightness variations induced by these processes is a 


powerful probe to understanding the underlying physical nature of sources. 

Systematic long term variability monitoring has been conducted at visible wavelengths 
since the beginning of the 20^ h Cen tury, espe cially following the discovery of the Cepheid 


Period-Luminosity relation (Leavitt fe Pickering! Il91 21. More recently, the advent of space- 
based infrared observatories extended these studies to the thermal infrared. At thermal 
infrared wavelengths emission is mainly due to warm astronomical dust heated by stellar 
sources or accretion processes. Infrared variability studies offer a unique window to char¬ 
acterize the temporal variations in stellar mass loss and accretion. Indeed, data from the 
InfraRed Astronomical Satellite (IRAS) and the Infrared Space Observatory (ISO), have 


van der Veen & Habine 

1990: 

Omont et al. 

1999: 

Whitelock et al. 

1991 

.), as well as the prop- 

erties of disk accretion 

m optically variable YSOs (e.g., 

Juhasz et al. 

2007; 

Abraham et al. 


2004 1. These pioneering studies were, however, limited to bright isolated objects in the 


Galaxy due to the low sensitivity and spatial resolution of the available thermal infrared 
instrumentation. 


The Spitzer Space Telescope (Werner et al. 


20041 ) and the Wide-held Infrared Survey 
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Explorer (WISE; Wright et ah 2010) for the Erst time provided the sensitivity required to 
repeatedly survey large areas of the sky to register the temporal variations of mid-IR sources. 
The observing strategy of the WISE mission, in particular, allowed the satellite to scan 
across each object in the sky at least 12 times every 90 minutes. The time series photometry 
obtained by this observatory are available as part of the recently released A11WISE catalog 
and have been u sed to characterize the light curves of short period variables sources in 


the Galaxy (e.g., [Klein et al. 2014). Spitzer, which operates as a facility telescope, allows 


monitoring with arbitrary timescales, constrained only by on-sky visibility as the observatory 
orbits around the Sun. Spitzer’s superior sensitivity and angular resolution (as good as ~2 
arcsec, a factor of three better than WISE) opened the possibility of monitoring the infrared 
variability of stars and disks in crowded areas like star forming regions in the plane of the 
Galaxy (e.g., Giannini et ah 120091 ). globular clusters in the Galactic halo, and the stellar 
population of the Galactic bulge (e.g., the ongoing Spitzer Carnegie RR Lyrae Program, 
monitoring stellar variability in 31 Galactic clusters and the bulge). Implementing Spitzer 
as a research telescope also rendered Local Group galaxies accessible for variability studies 
at infrared wavelengths. 

The _Spitzer Legacy Programs “Surveying the Agents of a Galaxy’s Evolution” (SAGE; 
Meixner et al.l 2006) and “Surveying the Agents of Galaxy Evolution in the Ti dally Stripped, 
Low Metallicity Small Magellanic Cloud” (SAGE-SMC; Gordon et ah 2011a), in particular, 
are ideally suited for studying the infrared variable population in the Magellanic Clouds, two 
of the closest Milky Way satellites. The scientific goals of these two surveys focus on the life 
cycle of baryonic matter as traced by dust emission, from its start in the Interstellar Medium 
(ISM) where new stars form, to the chemical enrichment of the ISM by the dusty winds of 
evolved stars. The two pr ograms imaged the Magellanic Clouds with both the InfraRed 
Array Camera (IRAC; [Faz io et al.l l2004h 3.6, 4.5, 5.8 , and 8.0 pm bands, and the Multiband 
Imaging Photometer for Spitzer (MIPS; R ieke et ah 2004) 24, 70, and 160 pm bands. Each 
galaxy was mapped twice in those surveys, allowing the detection of infrared variability over 
timescales roughly commensurable to the separation (several months to three years) between 
th e two ep ochs. An analysis of the SAGE Large Magellanic Cloud (LMC) data was presented 


m 


Vijh et ah 020091 ). Using the IRAC 3.6-8.0 /im and MIPS 24 pm data they found ~2000 


variable sources (out of over 4 million point sources in the SAGE catalog), the majority of 
which (~87%) were identified as evolved stars in the Asymptotic Giant Branch (AGB) phase. 
A small population of variable YSO candidates was also detected by Vijh et al. (i2009l ). 


In thi s work , we focus on the SMC, following a sim ilar approach to that developed by 


Viih et al.l ( 20091) . The SMC, with a distance of ~61 kpc ( Haschke et al.ll2012l: iHilditch et ah 


is characterized bv an even lower metallicitv than the LMC (~l/5 solar vs. ~l/2 

Russell & Dopita 

1992; 

Rolleston et al. 

1999, 

2003; 

Lee et al. 

2005)), below the thresh- 
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old, where the properti es of the ISM are expected to drastically change, as traced by th e 
low dust mass fraction fjEngelbracht et al.ll2005t iDraine et al.ll2007t ISandstrom et al.ll2010l) . 
The morphology of the SMC is characterized by a pronounced “Bar” with a “Wing”, and 
an extended “Tail” pointing toward the LMC (see Figure [T]). The latter structure is part of 
what is generally know n as t h e Mag ellanic Bridge, the H I filament connecting SMC with the 
LMC 


,e.g. 


McGee fe Newton 19861 ). The Magellanic Bridge is a clear manifestation of the 
tidal interaction between the two gal axies and of a possible close encounter ~200 Myr ago 
( Zaritskv fe Harris 2004 ; Harris 2007 ). The SAGE-SMC observations covered the full SMC 
including the Bar, Wing, and the Tail. We have supplemen ted this d ata with the pathfinder 
“Spitzer Survey of the Small Magellanic Cloud” (S 3 MC; Bolatto et ah 2007 ). to provide a 
third epoch photometry for the portion of the SMC Bar and Wing. While the SMC has 
been extensively monitored for optical variability as part of the MAssive Compact Halo Ob¬ 
jects Projec t (MACHO; lAlcock et al.lI1997T) and Optical Gravitational Lensing Experiment 


(OGLE-III; Udalsk i et al. 


microlensing experiments, this work is the first extensive 


survey of the SMC searching for variability at thermal infrared wavelengths. 


This paper is organized as follows. In Section [21 we describe the Spitzer SAGE-SMC and 
S 3 MC point source catalogs, as well as other optical and infrared data we used to identify and 
characterize the SMC variables. Section [3] outlines our methods for selecting and validating 
candidate variable sources. In Section [4] we compare the list of SMC variables to SMC 
sources classified with other methods, available in the literature. The demographics of the 
variables identified in this study is then discussed in Section 0 and a summary of our Endings 
is given in Section [HI 


2. The Data 


We utilized the data taken in three separate epochs with Spitzer during the cryogenic 
mission using the IRAC 3.6, 4.5, 5.8, and 8.0 pm bands, and the MIPS 24 pm band. As 
mentioned before, the data were collected as part of the SAGE-SMC and S 3 MC surveys: 


1. The Surveying the Agents of Galaxy Evolution in the Tidally Stripped, Low Metallicity 
Small Magellanic Cloud (SAGE-SMC; Gordon et al. 2 011a ) mapped the SMC in all IRAC 
(3.6, 4.5, 5.8, and 8.0 pm) and MIPS (24, 70, and 160 pm) bands in two separate epochs, 
which we identify as Epoch 1 and Epoch 2 for the remainder of this paper. These two 
epochs are separated by ~3 months (IRAC, data obtained in June and September 2008) 
and ~9 months (MIPS, acquired in September 2007 and June 2008), respectively. The 
observations covered the full SMC, including the Bar, Wing, and part of the Magellanic 
Bridge (the Tail; see Figured]), for a total coverage of ~30 deg 2 . 
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2. The Spitzer Survey of the Small Magellanic Cloud (S 3 MC; Bolatto et ah 20071) only 
mapped the Bar and Wing of the SMC (~8 times smaller coverage than SAGE-SMC, 
for a total of ~4 deg 2 ; see white contour in Figured]) in one single epoch (November 2004 
for MIPS and May 2005 for IRAC), preceding the SAGE-SMC survey by about 3 years. 
The S 3 MC IRAC (3.6-8.0 pm) and MIPS (70-160 pm) data were reprocessed using the 
SAGE-SMC pipeline as Epoch 0 data, providing a consistently reduced data set for all 
three epochs. In the following we refer to the reprocessed S 3 MC data as Epoch 0. 


The SAGE-SMC survey provides two catalogs containing the sou rces’ p hotometry for 
individual epochs, i.e. the IRAC Catalog and the IRAC Archive (see Gordon et al. 12011bl ). 
For onr study, we selected the more reliable Catalog over the Archive, favoring detection reli¬ 
ability and photometric accuracy over completeness. Similarly, we adopted the more restrie 


tive_MIPS 24 pm Catalog over the more complete, but less reliable Full List (IGordon et al. 
2011b ). The IRAC Catalog contains ~217,000, ~1.23 million, and ~1.13 million sources in 


Epoch 0, Epoch 1, and Epoch 2, respectively. The MIPS 24 pm Catalog includes ~5,600 
Epoch 0, and ~16,500 Epoch 1 and 2 sources. We did not use photometry from the MIPS 70 
and 160 pm bands due to the lower spatial resolution (~18" and ~40", respectively) relative 
to the shorter wavelength data. The IRAC and MIPS 24 pm data have angular resolution 
of 177-2" and 6", respectively. The Epoch 0 IRAC photometry has a limiting sensitivity of 
0.045, 0.028, 0.12, and 0.10 mJy in the IRAC 3.6, 4.5, 5.8, and 8.0 pm bands, respectively. 
The sensitivity of the Epoch 1 and 2 is about a factor of two lower due to the shorter in¬ 
tegration time. The MIPS 24 pm sensitivity is instead similar in the three epochs, and is 
~0.7 mJy. 

In addition to the source catalogs, our study also uses the SAGE-SMC Epoch 0-2 IRAC 
and MIPS 24 pm images for source quality validation. 


The MIPS data were consistently taken before the IRAC data with a period of no less 
than three months between observations. The MIPS and IRAC data were taken at different 
times within each epoch, which can allow us to use the data from MIPS to confirm a source 
detected as variable based on the IRAC data on a different temporal cadence. Table [Q 
provides information about when the data were taken by IRAC and MIPS for each epoch. 
The longest time between observations with IRAC and MIPS is nine months in Epoch 1, 
contrasted to the shortest time of three months in Epoch 2. Epoch 0 IRAC and MIPS data 
were taken six months apart from each other, but the time gap between Epochs 0 and 1 is 
nearly three years. 


The time interval between the different measurements is the main limiting factor to the 
typical period of the sources that we are able to detect in our dataset. The detectability 
of a source as a variable will depend on how well the sampling of the light curve matches 
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the period of variability: our study is most sensitive to sources which our instruments are 
able to detect in different phases of their variability period when we collect our data. These 
types of sources include long period variables that have not completed a full phase of their 
variability period during our study. Sources that are in the same phase of their variability 
period when our data were collected, will not be detected as variable. As a consequence, 
our study is biased toward redder sources with a long period of variability (of the orders of 
months to years) with high amplitudes of variations in the thermal infrared. 


2.1. Ancillary Data: Optical and Near-IR Point Source Catalogs 


The near-IR JHK S data from the combi ned Two Micron All-Sky (2M ASS; ISkrutskie et al. 


2006 1 and the 2MASS 6X Deep (6X2MASS; Cutri fc 2MASS Team 2004h Point Source Cat 


alogs are the integral part of the SAGE-SMC IRAC catalogs (see [Go rd on et al. 20 11b). The 
SAGE-SMC catalog is cross-matched in the SAGE-SMC database with th e opt ica 


the Magellanic Clouds Photometric Survey (MCPS, UB VI bands: IZar i 


.skv et ah 


data from 
2002) and 


Udalski et al. 


2008). 


the Optical Gravitational Tensing Experiment (OGLE-111, VI bands; 

We considered IRAC and optical sources a match when a distance between their positions 
was <1". 


The ancillary data were used to construct sources’ spectral energy distributions for a 
wide wavelength range, from optical to 24 /irn. We adopt the V- and /-band photometry 
from OGLE-III when available, from MCPS otherwise. 


Variable Star Identification 


We identified the SMC variable sources following the procedure described in Vijh et al. 


(120091 ) for the LMC. We first matched the individual SAGE-SMC Epoch 0, Epoch 1, and 
Epoch 2 catalogs to extract a time series for each source detected in all epochs. We then 
applied variability criteria based on the flux change between epochs to select sources with 
measurable change above the photometric sensitivity. We finally performed quality checks 
aimed to remove spurious detections in order to have a reliable catalog of SMC infrared 
variables. 
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3.1. Interepoch Matching 

In the initial step of our variable stars selection process, we cross-matched single epoch 
IRAC and MIPS 24 //rn catalogs. Accurate matching between epochs is crucial to ensure 
that we are measuring changes in flux for the same sources across multiple epochs. The 
match was performed using the CasJobs interface within the Mikulski Archive for Space 
Telescopes (MAST). CasJobs is based on a SQL database. 

We started by generating a list of all the sources that had coordinates within 1" of 
each other within at least two of the three epochs of the SAGE-SMC IRAC data. We then 
matched the resulting IRAC lists with their corresponding MIPS 24 /jrn data for each epoch, 
again requiring a match within 1". The requirement to identify common sources across 
multiple epochs of data significantly reduced the number of sources available for our study. 
For instance, the largest of our matched lists, containing sources with matches in Epoch 1 and 
2 (the largest spatial coverage of the observations), contains just over 800,000 sources with 
valid 3.6 /jrn photometry, compared to over 2 million sources in each single epoch (Epoch 1 
and Epoch 2) catalog. The source counts for other IRAC and MIPS 24 /jrn bands and pairs 
of epochs can be found in Table [2j Note that our method of searching for variables would 
miss novae and supernovae if the progenitor star was not detectable in our survey before the 
outburst. 


3.2. Variability Criteria 


We identified candidate variables out_of our inter-epoch matched catalogs on the basis 
of the variability index V defined by Viih et al.J (120091 ) as the error-weighted flux difference 
between each pair of epochs: 


V = 4 = (1) 

s] a l + T 

where _/) and fj are fluxes in a certain band for a single source in epochs i and j, respectively; 
and (Jf z and aj j are respective flux uncertainties. This index was calculated for each possible 
pair of the three available epochs ( i,j = 0,1,2 and i ^ j ) in each band, so that a single 
source can have up to 15 variability indices if it has valid fluxes in all the IRAC and MIPS 
24 /jrn bands across all three epochs of data. Note that a positive V would correspond to 
an increase in flux between the two epochs, while a negative index would indicate that the 
source decreased in flux from one epoch to the other. 







We considered a source to be variable if \V\ >3 (equivalent to measuring flux changes 


with at least 3a sensitivity) in at least two neighboring bands, again following Vijh et ah 

fhoOhlh 


I^3.6ij| 

> 3 

and 

imq 

> 3 

IvLbJ 

> 3 

and 

U.Sij 1 

> 3 


> 3 

and 

U.Oij 1 

> 3 

1 ^ 8.0 J 

> 3 

and 

\Vl4 ts \ 

> 3 


with i 7 ^ j and i,j = 0 , 1 , 2 


( 2 ) 


where i and j are two of the three available epochs. We also required that the variability 
index among neighboring IRAC bands have the same sign, since all the IRAC bands were 
observed at essentially the same time within each epoch and we expect the change from 
one epoch to another to be similar for each band. Failing this condition is an indication 
that the flux in different bands may come from mismatched sources, sources whose flux is 
contaminated by a nearby or blended source, or by the presence of an artifact in our data. 
The exception to this rule is for the 24 /mi photometry, for which we do not require the 
same sign variability; the MIPS data were taken at different times than the IRAC data and 
correspond to a different variability phase. We finally excluded all sources with magnitudes 
larger than 15 at 3.6 /mi and 4.5 /mi, 13 at 5.8 /mi, and 12 in the 8.0 /mi IRAC bands; 
the photometric uncertainties increase for fainter sources, making flux measurements less 
reliable f Gordon et al. 2 011 b ). 


Note that to be considered as variable star candidates, sources do not need to be variable 
between all three pairs of epochs (only ~ 6 % of all candidate variable sources have |R| > 3 
in all three epoch pairs). Sources do not need to have valid photometry in necessarily all 
IRAC and MIPS bands as long as the sources have enough valid fluxes to demonstrate their 
variability. 


The results of our search for variable stars is summarized in Table [2] and the distributions 
of the variability indices for all bands and epoch pairs are shown in Figured The variability 
indices for IRAC bands between Epochs 0 and 1 and Epochs 0 and 2 have a larger width (as 
shown in Figure [2D than the variability indices between Epochs 1 and 2 . This is a consequence 
of the longer integration time used in the S 3 MC survey (reprocessed into SAGE-SMC Epoch 
0 ), leading to smaller flux uncertainty and hence a better sensitivity for flux changes measured 
by the variability index. Conversely, the uncertainly in the MIPS 24 /mi photometry (and 
hence the spread in the variability index distribution) was similar between the SAGE-SMC 
and S 3 MC datasets. We found 825 sources that meet our criteria for variability in at least 
two consecutive bands and pairs of epochs. 









9 


3.3. Source Inspection 

We have visually inspected each variable star candidate in all available bands and epochs 
using thumbnail images extracted from the SAGE-SMC IRAC 3.6, 4.5, 5.8, and 8.0 /mi, and 
MIPS 24 /mi band mosaics. The purpose of the visual inspection was to remove candidates 
whose photometric reliability was compromised by poor pixel quality, optical or electronic 
artifacts, un-rejected transients such as cosmic rays, or contamination from a neighboring 
source. At the end of this process we were left with 814 candidates that meet all our criteria 
for variability, photometry, and image quality. 


3.4. The SMC Infrared Variables Catalog 

Our final catalog of infrared variables in the SMC is available in full within the online 
version of this paper. The catalog lists the basic astrometric and photometric properties 
of each variable candidate, including IRAC and MIPS designations, fluxes and magnitudes 
with uncertainties in the optical (MCPS and OGLE-III), near-IR (2MASS/6X2MASS), and 
rnid-IR (Spitzer IRAC and MIPS 24 /mi) bands, as well as the variability indices V for 
each pair of epochs. The catalog also includes possible classification of each source based 
on matches with the literature (see Section [4] below). The complete list of columns in the 
catalog of SMC variables is provided in Table [3] 


4. Source Classification 

In order to understand the demographics of 814 newly-identified variable star candi¬ 
dates, we have cross-matched our source list with a number of published catalogs of SMC 
sources. These catalogs were constructed using different methods and include sources classi¬ 
fied with varying levels of confidence. We required a distance between a variable source and 
its matching source from literature to be <1". A detailed discussion of the catalogs we used 
to determine the nature of variable sources and the matching results follows. The summary 
is provided in Table [4] 

For several SAGE-SMC variable sources we found multiple matches in literature. Re¬ 
solving ambiguity in source classification is out of scope of this paper, thus we provide all 
available classifications in our catalog (columns 19 to 21 in Table |3|). 
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4.1. Evolved Stars 


Over 80% of the LMC infrared variables identified by Vijh et ah (12009! ) are evolved stars: 
a similar result is to be expected in our SMC search. Stars of all masses, once they leave 
the main sequence, enter one or more phases when their convective envelopes expand by 
several orders of magnitude, while their effective temperatures become as low as ~3,000 K. 
In those conditions, these stars cross the Long Period Variables (LPVs) instability strip in 
the Hertzsprung-Russell diagram, where radial pulsations are modulated by periodic changes 
in the opacity and ionization of internal stellar layers. These pulsations are observed as lu¬ 
minosity variations with period ranging from several tens to thousands of days. According 
to the shapes, amplitudes, and periods of their optical light curves, LPVs are classified as 
Mira, Semiregular, or Irregular variables (jSamus et al. 2012|), with Miras typically having 
larger amplitude luminosity variations. LPVs are also characterized by intense mass loss, 
which leads to the formation of dusty circumstellar envelopes. The circumstellar envelopes in 
turn give rise to strong infrared excesses above the photospheric continuum. While the vari¬ 
ability properties of LPVs are well characterized in the visible^ their mid-infr ared variability 
has only recently been studied in a systematic way (e.g., R iebe l et ah 2010 and references 
therein). 

Low and intermediate mass stars (M < 8 M 0 ) will become LPVs twice in their life: as 
stars in the Red Giant Branch (RGB, characterized by H nuclear burning in a shell around 
an inert He core), and in the Asymptotic Giant Branch (AGB, characterized by alternate 
H and He burning in a shell around an inert C/O core). While the AGB phase is shorter 
than the RGB phase, AGB stars tend to be more luminous, show larger variability, and 
have more intense mass loss. As a consequence, our search is biased toward discovering 
AGB rather than RGB variables. AGB stars are subjected to third dredge-up events, in 
which elements synthesized during He-shcll burning events are transported toward the stellar 
surface. The end product is the injection of trace elements (the so-called s-elements) in the 
stellar atmosphere, uniquely produced during these He-shcll burning events. Third dredge- 
up events ultimately lead to a gradual increase in the photospheric C/O abundance ratio 
that, for a star with a main sequence mass in the ~ 1.5 - 4 M 0 range (for solar metallicity; 
Straniero et al.lIl995l . 1997}), leads to the formation of a so-called carbon star. Carbon stars 
are characterized by a C-rich atmospheric and circumstellar dust and gas chemistry. Stars 
below this mass limit never experience enough third dredge-up episodes to become C-rich 
AGBs, while stars above the upper limit destroy the newly generated C in a proce ss called Hot 
Bottom Burning (HBB; ISmith fe Lambertlll985l; iBoothrovd fc Sackm an nil 19921 1 . Stars that 
undergo HBB processes remain O-rich AGBs and tend to have much higher luminosities than 
low mass O-rich AGBs. Note that in the low metallicity environments, like in the Magellanic 
Clouds, stars have a lower initial O atmosphere abundance, and a smaller number of dredge- 
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up events are required to push the C/0 ratio above unity. In this case the lower mass limit 
for the product of C-rich AGBs is somewhat reduced and a larg er population of C-rich AGBs 
is to be expected (IFerrarotti fe Gailll2006t lYentura et al.l 120121) . 


The mass loss rate dramatically increases toward the end of the AGB phase (jWillson 


20001); in this late stage, the circumstellar envelope of AGB stars becomes optically thick 
and completely obscures the star in optical wavelengths. Stars of this kind are often referred 
to as extreme AGBs and can be either O-rich (in which case they may show OH/IR maser 
activity) or C-rich. A handful of the most extreme A GB stars are responsible for the large r 


fraction of dust injection into the LMC and SMC ISM (Boyer et al. 2012yRie 


which these dust contributors tending to be prevalently C-rich (iBover et al 


pel et al.ll2012f) . 

2 01lh . 


High mass stars, on the other hand, cross the LPV instability strip as Red SuperGiant 
(RSG) stars. RSG stars are characterized by a higher luminosity than low mass AGB stars 
(but comparable to intermediate mass HBB O-rich AGB stars) and a similar mass loss rate. 
RSG stars tend to have much smaller populations than AGB and RGB stars due to the 
Initial Mass Function being biased toward lower masses, and the shorter life span of RSG 
stars. As such, RSG stars are not large contributors to the ISM and we should not expect 
to find many of them among our candidate variables. 


We matched our SMC variable sources with the color-selected evolved stars fromlBover et al. 


(2011) and with the LPVs identified in the OGLE-III project by Soszvhski et ah (2011) based 
on the characteristics of their light cu rves (and then classified as O-rich or C-rich on the basis 
of their optical and near-IR colors). Boyer et al. 1201111 differentiate the stars according to 
the evolutionary phases (RGBs, O-rich AGBs, C-rich AGBs, extreme AGBs, and RSGs), 
while ISoszynsk i et al. 120Ilf) only differentiate between RGBs, O-rich AGBs, and C-rich 
AGBs (along with providing details on the variability type). We also matched our sample 
of SAGE-SMC variables with the catalog of massive stars from Bonanos et a l. 02010) that 
includes spectroscopically identified RSG stars. As expected, most of the matches with our 
variable candidates are AGB stars, of which 82 are O-rich and 497 are C-rich AGBs. Of 


those O- and C-rich AGBs, 251 are classified as extreme AGBs according to IBover et a l. 


(120111) . the vast majority o f which (85%) are C-ric h according to 


S oszvh ski et ah (1201 ll) . 


For the non-extreme AGBs, Boyer et ah (2011) and ISoszynski et ah (20111) are generally in 
agreement, with ~16% mismatch within the O-rich category and less than 3% disagreement 
for the carbon s tars. We only found 6 RGB and 23 RSG stars in our study, but it is worth 
noting that the Boyer et ah- (2011) classification scheme cannot distinguish between RGB 
stars and AGB stars with luminosity below the tip of the RGB. 


Boyer et ah (2011) has also three categories that we have not considered in our analysis. 


Sources labeled as “Far InfraRed objects” (FIR) are characterized by increasing flux at 8 








































































12 


and 24 pm, even larger than what is expected for extreme AGB stars. The fact that the 
FIR sources are mostly found in proximity of star forming region supports the hypothesis 
that they may be misclassified YSOs rather than evolved stars, but without spectroscopic 
analysis it is not possible to rule out the possibility that our FIR sources are evolved stars 
with unusually thick dusty envelopes. For this reason we have left FIR sources as unclassified 
in our survey. The “anomalous” aO-ri ch and aC-r ich AGBs are O-rich and C-rich star, 
respectively, with peculiar colors in the Boyer et ah 1201111 analysis. Rather than keeping 
these intermediate c lasses, we h ave relied on the chemical classification provided by the 
OGLE-III analysis in Soszynski et al.J 12 01111 . which is based on a combination of optical and 
near-IR colors. We found that most (18 out of 27) aO-rich AGBs are confidently described 
as O-rich AGBs in the OGLE-III catalog, with the remaining being classified as C-rich AGB 
stars, while the two aC-rich AGBs are evenly split between C- and O-rich AGBs in OGLE-III. 

Due to their lower brightness, we were unable to calculate the V 24 index of RGB stars 
for any of our intervals of variability. RGB stars also have |Vs.o| < 3 across all variability 
intervals. One source, SSTISAGE1C J004906.89-724309.7, has contrasting classifications 
among different catalogs. Soszynski et ah (2011) lists the source as a RGB star with O-rich 
chemistry and Bo yer et ah| (2 0111) classifies the source as a C-rich AGB star. This source is 
the only RGB candidate we found that varies in the 5.8 pm band (all other RGB stars show 
variability in the 3.6 and 4.5 pm IRAC bands and have |V 58 1 < 3 within all epochs). 

Two RSGs have c ontrasting mat c hes in literature. SSTISAGE1C J010302.41-7 2 0152. 9 
is classified as RSG by Bonanos et ah ( 2Qloh and as O-r ich AGB star by Bover et ah ( 2011 ). 
Given the higher reliability of the Bonanos et ah (120101 ) classification that used full spectra 
rather than broadband colors, SSTISAGE1C J010302.41-720152.9 is likely to be an under- 
luminous RSG star. This source meets our variability criteria primarily in the IRAC bands 
between Epochs 0 and 1, but the MIPS 24 pm band has the largest variability of all RSG 
sources with V 24 = —15.80 ( betwee n Epochs 0 and 2) and —12.92 (between Epochs 1 and 
2 ). Note, however, that ISoszvhski et ah (20111 ) is not attempting to differentiate AGB from 
RSG stars. RSG stars appear in the OGLE-III catalog as Semiregular LPVs with O-rich 
chemistry. This_suggests that this source may be a RSG based on the K s luminosity cut in 
Bover et ahl (2011J). Interestingly, this source is the most variable source between Epochs 1 
and 2 in terms of consecutive neighboring bands with \V\ > 3. 

Finally, we have identified a candidate planetary nebula (PN), matched from a catalog 
kindly provided by G. Jacoby (2014, private communication). However, this sam e sou rce 


(SST ISAGE1C J005530.37-725021.9 ) is also cross-m a tched with a YSO from iBolatto et al. 


(2007) and a background AGN from Kozlowski et ah (20111). The variability of this source is 
detected only between Epochs 0 and 1 and Epochs 0 and 2 in the 5.8 and 8.0 pm bands. This 
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suggests that this object may be changing brightness on very long time scales (comparable 
to the ~3 years separation between the S 3 MC and SAGE-SMC surveys), pointing toward 
changes in accretion rate as a possible mechanism for the infrared variability. 


4.2. Cepheids 


Intermediate mass stars (4-8 M 0 ) cross the Cepheid instability strip during the He-core 
burning stage, while they describe a pronounced “blue loop” on the Hertzsprung-Russell 
diagram. Cepheids are radial pulsators similar to LPVs, but their typical variability period 
is between 1 and 80 days. Ceph eids have F and G spectral types; typically they do not 


have substantial mass loss (but see lMarengo et al.ll2010l;lBarmbv et al.ll201ll ; Matthews et al. 


2012 ) and their intrinsic brightness makes them easily detectable in the SAGE-SMC survey. 
Given their imp ortance as d i stan ce indicators, thanks to their period-luminosity relation 
(Leavitt law, Leavitt & Pickering [1912), Cepheid variables have been extensively studied 
not just in the M ilky Way, but also in local group galaxies (includ ing the SMC) at visible 


2010j; Scowcroft et al 


2011 


2013) 


wave l engths (e.g.. lFreedma.ii et al.ll200ll) . as well as in the infrared fjNgeow fe Kanburl 12008 


To search for Cepheids, we matched our catalog of candidate variables with the source 
lists in the OGLE-111 database (where, again, the classification is based o n the ch aracte ristics 
of the light curves). We matched both C lassic al Cepheids flSoszvhski et ahl 2010a') and 
the lower mass Type 11 Cepheids (iSosz ynski et al. l2010bf l. We found 70 Cepheids in our 
catalog (9% of our sample), with 67 of sources being classified as Classical Cepheids, and the 
remaining 3 being classified as Type If Cepheids. Our Cepheids show the highest variability 
index in the 3.6 and 4.5 /jrn bands, but generally have |C| < 5 across all intervals of variability 
in all bands. 

We also matched our catalog with the list of OGLE-111 RR Lyrae in the SMC from 
Soszvriski et al. (2010q). We did not find any matches, most likely due to the lower ampli¬ 


tudes and luminosities that RR Lyrae variables have in the thermal infrared (Madore et al. 
2013h . 


4.3. Early Type Stars 

Early type stars (high mass stars with typically O and B spectral types) are characterized 
by several different kinds of variability, in some cases due to periodic pulsations (e.g., radially 
pulsating f3 Cephci and non-radially pulsating a Cygni stars), in others due to episodic 
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outbursts (e.g., Luminous Blue Variables, whose prototype is /3 Dor in the LMC). In other 
cases, variability in early type stars is caused by accretion in binary systems and can manifest 
itself at wavelengths ranging from the X-rays to the infrared (e.g., X-ray binaries). Early 
type stars tend to be rare, even more so the ones that show detectable variability. 


To search for early type variables in the SMC, we have cross matched our varia bl e star 


candid ates list with the catalog of spectroscop i cally confirmed massive stars from lBonanos et al. 


(1 2010 ) and dusty OB stars from [S heets et ah (1201311 . We found a total of 23 sources (~3% of 


—^ / 

Sheets et al. ( 

2013 

- - ^ - --- - 

). Four of the early spectral type stars identified by 

Bonanos et al. 

(2010) 


are also matched with the Sturm et al. (20131) catalog of SMC X-ray sources. Those sources 
are likely accreting X-ray binaries. 


4.4. Young Stellar Objects 


According to the Hartmann & Kenyon (19851) model of accretion outbursts in FU Orio- 
nis objects, infrared variability in YSOs is generally associated with changes in the accretion 


rate modulated by the protostellar disk rather than precessing jets (Herbig et al 
which are rather invoked to explain optical variability. 


2003), 


Observations at different wavelength regimes allow a detection of YSOs at different 
stages of their evolution. The SED of a YSO changes when it evolves. Stage 0 protostars 
are in the early protostellar collapse stage with the majority of the source mass residing 
in the infalling envelope. The SED is dominated by the far-IR, because the protostar and 
possibly a disk are deeply embedded in a circumstellar envelope that absorbs the forming 
star’s photospheric emission, warms the dust and re-emits it in the far-IR. Stage I is a later 
stage of protostellar collapse. Stage I YSOs display very broad SEDs that peak near 100 
/im; they emit the bulk of their radiation at the mid-IR to far-IR wavelengths. The envelope 
mass is similar to the mass of the central protostellar object. Stage I YSOs have well- 
developed accretion disks and their envelopes have bipolar cavities excavated by outflows. 
Stage II YSOs are characterized by the presence of optically thick disks, bipolar outflows, 
and possibly the remains of a tenuous infalling envelope. Stage II YSOs emit mostly in 
the near-IR to mid-IR. Stage III YSOs, which are mostly detected at optical to near-IR 
wavelengths, have optically thin disks or no disks, and may still have optical bipolar jets. 


The mid-IR Spitzer data are sensitive to Stc 

ige I and Stage II YSOs. Prior to Spitzer , 

there was only one YSO identified in the SMC 

f Gat lev et al. 

1982 

). In the galaxy-wide 

study based on the Spitzer SAGE-SMC data, 8 

ewilo et al. (2 

013) identified ~1000 YSO 
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candidates which are categorized in to “possib le YSO candidates” and “high-reliability YSO 
candidates.” An earlier study by Bolatto et al. 020071 ) that used the S 3 MC data (and thus 
concentrated on the main body of the SMC) found ~300 YSO candidates; this list is less 
vetted than t hat o f lSewilo et al.l (120131) and potentially more contaminated by non-YSOs. 


Chen et ah (20141 used the SAGE-SMC data to conduct a photometric search for YSOs in 
the SMC Tail and identified 26 massive YSO candidates. More detailed studies on YSOs 


in the SMC co ncentrated on individual star-forming regions: NGC346 (ISimon et al.l 1200711 


and NGC602 ( Carlson et al. 1120071: ICarlson et al.l 12011 


Gouliermis et ah 


2007), identifying 


~150 additional YSO candi dates; their YSO lists include low luminosity sources missed by 
galaxy-wide studies. Sewilo et ah (120131) used the relatively conservative brightness cuts to 
remove typically dimmer background galaxies from their YSO list, and thus selecting the 
brightest, most massive YSOs. The detailed inspection of source morphology that would 
help to distinguish between background galaxies and YSOs can only be done on smaller 
scales. To date, Spitzer studies identified over 1,200 individual YSO candid ates in the SMC; 
33 sources are conf i rmed spectroscopically as bona fide YSOs (Ivan Loon et al.l 120081. l2010l; 
Oliveira et al. 2011 . 2013 b 


Operating at longer wavelengths, the Herschel Space Observatory (IPilbratt et ah 20101 


is sensitive to lower dust temperatures than Spitzer , allowing a detection of YSOs at the ear¬ 
liest evolutionary s tages that had been missed by Spitzer. Herschel is sensitive to Stage 0 
and Stage I YSOs. ISeale et al.l (1201411 used the data from the Hersc hel Key Project HER¬ 


ITAGE (“Herschel Inventory of the Agents of Galaxy Evolution”; Meixner et ah 120131 ). 


which imaged the Magellanic System at 100, 160, 250, 350, and 5 00 pm, to identify pop¬ 
ulations of sources in the LMC and SMC. In the SMC, Seale et ah ( 2014 1 identified ~700 
YSO candidates, which are categorized into “probable YSO candidates” and “possible YSO 
candidates”; 24% of these are newly-identified YSO candidates. 

We crossmatched our list of candidate SMC variables with the YSO candidate lists 
from literatu re to identify i nfrare d vari able YSOs. We fo und 4 4 matches in 


(201311. 13 in 


Carlson et ah 


Bolatto et al.l (1200711 . 10 in ISea 


(20111). and one match in each 


e et ah (2014b 


van Loon et ah 


m 


Oliveira et ah 


Sewilo et ah 


20081). Ivan Loon et al.l (2010), 


(2013), 3 in 


and Oliveira et ah (20 1111 (all 3 refer to the same source). Many sources appear on more 


than one YSO list from literature; our list of candidate SMC variables includes 47 individual 
YSO candidates. 

However, out of 47 sources matched to YSO candidates, 11 were also matched to other 
classes of sources. Five sources were matched to extreme AGB stars (IBover et al.l 201111. 
tending to have very high variability indices at the 24 pm band (| V 24 1 > 20) b etwee n Epochs 
1 and 2. The remaining sources were matched to the early type stars from Bonanos et ah 
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(2010I), and the multiple-match PN described in Section 14.11 above. 


Of the 984 YSOs in Sewilo et al. 020131 ). we match 44, suggesting a lower limit of ~4 
SMC YSOs th at a re variable in the infrared with our sensitivity and timescale coverage. 


Sewilo et al.l (1201311 estimated physical parameters of 22 of th ese sources based on th e results 


of the SED fitting with the YSO radiative transfer models of (Robit aille et al. 2006). All 22 
sources are classified as Stage 1 YSOs, indicating that the infrared variability may be tracing 
the accretion variability. 


4.5. Background Galaxies 


The OGLE-111 survey has been used to search for optical variability of AGNs and 
quasars (QSOs) as part of the Magellanic Quasars Survey (IKozlowski et al.l 201 ll. 20131). 
These sources are clearly background objects and not part of the SMC stellar population. 
Still, Ending background galaxies within our survey is a useful complement to AGN and 
QSO optical studies, which typ ically sho w flux changes between ~3-30% in the visible, over a 
period of 10-100 days (jKellv et ah 2009). The detection of variability at Spitzer wavelengths 
would indicate the presence of a non-thermal component in the characteristic active galaxies ’ 


mid-IR spectral “bump” (e.g.. Ilmpev &; Neugebauerlll988l: iNeugerbauer &; Matthewslll999 
and references therein). Furthermore, identification of distant quasars projected within the 
held of view of the SMC can be a useful tool to determine the proper motion of this galaxy. 


We matched our catalog with IKozlowsk i et ah (201 ll. 12013 ) lists of AGNs and QSOs, 
which have been identified through spectroscopy and light curve morphology collected over 
a time span of 12 years. We found 8 sources ma tching the IKozlowski et al.l (1201 ll) AGN 
list, and one matching the Kozlow ski et ahl (1201 34 QSO sample. The latter source is also 
matched to an X-ray source from the Sturm et ah ( 2013 ) catalog (SST1SAGE1C J011139.57- 
725031.6). As mentioned above, one of the variable sources matched to the AGN (SST1S- 
AGE1C J005530.37-725021.9) is also matched to a PN (G. Jacoby 2014, private communi¬ 
cation), and a YSO (Bolatto et al. 2007). 


4.6. Other Catalogs 


Finally, we have also matched our sample of variable candidates with other source lists 
with peculiar spectral properties , namely WhOOO X -ray sources from lSturm et ahl (120131) and 
a catalog of radio sources from Wong et ah J2012 1 (457 sources detected at 3 cm and 601 
sources at 6 cm). We found a match for 6 X-ray sources, the majority of which are associated 
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with early type stars (hence most likely X-rav binaries), and one (SSTISAGE1C J011139.57- 
725031.6) in common with iKozlowski et al.l (120131) QSOs. We o nly identified one 6 cm radio 
source as an infrared variable (less than 1% of I Wong et ah, 2012 catalog); this variable Spitzer 
source (SST I SAGE IC J005212.88-7308 5 2.8) i s also matched to the evolved stars from the 


Bo ye r et ah (201 ill and Soszvhski et ah (2011) catalogs (see discussion in Section I5T3T) . 


4.7. Unclassified Objects 

At the end of our classification process, we remained with 59 sources (~7% of the total) 
that are not matched to any of the catalogs described above. These unclassified sources tend 
to be fainter (by ~1 magnitude on average) than the classified variables and lack 24 /irn 
photometry (less than 20% of the unclassified sources have 24 /irn photometry in at least 
two epochs). The unclassified sources that do have at least one epoch detection at 24 /irn in 
most cases possess an SED resembling the one of YSOs (see discussion in Section IQ]) . 

Many of the unclassified sources are located in the SMC Tail (see Section [5]), as ex¬ 
pected since the Magellanic Bridge is not covered by many of the surveys we used for source 
classification. 


5. Discussion 

Figure El shows the distribution of the absolute values of variability indices in the 3.6 /irn 
band between Epochs 1 and 2, separately for each class of sources. We chose this wavelength 
and pair of epochs because this is the combination providing the largest number of sources 
for which a variability index is available (mainly due to the higher sensitivity of this band 
and full coverage of the SMC in Epochs 1 and 2). The histograms for other bands and 
combinations of epochs are similar. 

Figure [3] shows that most of the sources identified as variables according to the criteria 
delineated in Section I3T21 have a 3.6 /irn variability index below the minimum threshold of 
3. The fact that these sources are still classified as variables indicates that high variability 
is still detected in at least two other bands or in another pair of epochs. The exception are 
extreme AGBs with an average variability of 5.06, well above the threshold, and Cephcids, 
that have a median absolute variability index of 2.85, still significantly higher than other 
classes. It should be noted that these sources are the ones with the longest (extreme AGBs) 
and shortest (Cepheids) periods of variability among the type of variables we matched: this 
increases the chances of detection by our method by reducing the possibility of temporal 
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aliases with the inter-epoch matching in our datasets. 

A different way to look at the same behaviors is shown in Figure [U This figure plots the 
median value of the absolute variability index \V\ for each class, each band, and each pair 
of epochs, separately. The vertical bars represent the standard deviation in the distribution 
of \V\. The variability indices tend to be larger when calculated with respect to Epoch 0, 
in part because this epoch corresponds to the more sensitive observations obtained with the 
S 3 MC survey. The much shorter time interval between the two SAGE-SMC epochs (Epochs 
1 and 2), however, does play a role for some classes of sources. Cepheid stars, for example, 
show almost the same variability index between Epochs 1 and 2 as in the other two intervals, 
possibly because of their very short period which is equally sampled in all three intervals. 
RSGs, early type stars, YSOs, and galaxies, on the other hand, show very little variability 
between Epochs 1 and 2, even in bands (3.6 and 4.5 /ini in particular) where they have high 
variability in the other two intervals (Epochs 0 and 1 and Epochs 0 and 2). O-rich AGBs 
and YSOs have the largest spread in variability indices in the 24 /irnband; in all other bands 
the sources showing the largest standard deviation are instead the extreme AGB stars. 


Figure 0 shows the spatial distribution of the individual variables over the SAGE-SMC 
8.0 /irn map. The figure shows that most of the variables are concentrated in the SMC Bar 
and Wing, due to the larger stellar density in these parts of the galaxy. Young sources (YSOs 
and moderately massive objects such as early type stars, RSGs, and Cepheids) are over¬ 
represented in the Bar, which is where most star formation is observed to occur. However, 
one YSO is detected in the Tail (SSTISAGE1C J014915.82-743225.5), which indicates the 
presence of acti ve st ar formation sites in the bridge con necting the L MC w ith the SMC, as 
suggested by Harris (12007 ) and recently confirmed by Chen et ah (2014). Evolved stars 
are found at larger distances from the Bar, including the Tail, as expected for their lower 
mass and longer life span, which allows them to diffuse farther from the areas where they 
are born. Note a large number of symbols for extreme AGB stars overlapping with C-rich 
AGBs, which indicates that the extreme AGB stars in the SMC have in most cases a carbon- 
rich chemistry (at least in the areas of the Bar and Wing, where OGLE-III survey data are 
available). Finally, it should be noted that a large fraction of the variable sources in the 
SMC Tail are unclassified: this is a consequence of the lack of coverage in that area in the 
various source catalogs that we have used for source classification. 


5.1. Color-Magnitude Diagrams 

In Figures [6}i-d, we have plotted our candidate infrared variables in a number of color- 
magnitude diagrams (CMDs), which allow a better understanding of the nature of these 
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objects and can be used to search for anomalies (sources with an assigned class that show 
unique spectral properties for sources of that type). In all plots, the entire SAGE-SMC 
Epoch 1 catalog is shown in greyscale as a source density plot (a Hess diagram). 


Figure [0}i shows the near-IR CMD used for source selection by Boyer et al. 12 01 11 1. 
Indeed, in this diagram the evolved variable stars of different classes appear well sepa rate d, 
as a consequence of our source classification process heavily relying on thelBo ve r et al. 120111) 
catalog for evolved stars. In particular, the diagram demonstrates the neat color separa tion of 


extre me AGB stars from “regular” carbon stars (again a consequence of using the Bo yer et al. 


20111 catalog for source classif ication): note, however, how many of these sources are still 


classified as C-stars by Soszynski et ah (20111). Another prominent feature in this diagram 
is the vertical sequence of Cepheid stars (blue circles), which is well isolated from the redder 
evolved stars, and from the much bluer early type stars. The YSO sequence is also well 
separated and shows how sources with larger infrared excesses are also brighter in the K s 
band. Of particular interest among the outliers in this diagram is a very red (J — K s ~ 1.37) 
source (SSTISAGE1C J005107.19-734133.3) classified as a Type II Cepheid in the OGLE-III 
database. This source is discussed in detail in Section 15741 Figure [d}i also shows a variable 
source (SSTISAGE1C J010921.65-712435.0) along the YSO sequence that is matched to 
both the YSO and C-rich AGB catalogs. The OGLE-III light curve for this variable is 
consistent with being a Mira, which would support its classification as a carbon star, despite 
its infrared colors being more consistent with the source being a YSO. The anomalous colors 
of this source could be explained by a photometric blend of two different sources, or by 
the presence of a circumstcllar disk around a pulsating AGB star, which would mimic the 
spectral energy distribution of a YSO. Further analysis is required to confirm the nature of 
this source. 


A similar separation between different classes of SAGE-SMC infrared variables is shown 
in panels b and c of Figure |HJ The low and intermediate mass evolved star sequence is 
dominated by (in order of increasing brightness and redder colors) RGBs, O-rich AGBs, C- 
rich AGBs, and extreme AGBs (most of them classified as having a C-rich chemistry). The 
RSGs have infrared colors similar to O-rich AGBs (due to the similar dust composition), but 
generally higher luminosity. Note an isolated group of O-rich AGBs with [3.6] < 10, having 
a comparable brightness to RSGs, but moderately redder colors (by up to J — [3.6] ~ 1 mag, 
indicating higher mass loss rate): these are likely candidate intermediate mass AGB stars 
undergoing HBB. In these diagrams, the YSOs are better isolated from the main sequence 
and the lower giant branch, but the brightest members of this group (likely high mass YSOs) 
still have infrared colors and magnitudes analogous to the reddest extreme AGB stars. 

Figure [6ji is the CMD with the smallest number of sources due to the lower sensitivity of 
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the MIPS 24 //rri band. However, it provides a better separation of the YSOs from all other 
classes, including background galaxies (the latter being concentrated in the large clump with 
[8.0] > 12 and [8.0] — [24] > 2.5). Of all the infrared variable active galaxies identified in our 
search, two are detected at 24 /jm and are located at the upper end of this clump. 

We further investigate the infrared properties of candidate variables identified in our 
study by analyzing in detail their complete SEDs. This is discussed in the following sections. 


5.2. Spectral Energy Distribution 

Figures [7] to [13] show SEDs of a representative sample of variable sources of differ¬ 
ent types. The UB (MCPS), VI (OGLE-III if available, MOPS otherwise), and JHK S 
(2MASS/6X2MASS) data are plotted as black circles. All three-epoch IRAC and MIPS 
24 /jrn photometric data are plotted and color-coded as indicated in the plots. The data 
points from individual surveys (representing separate epochs of observations) are connected 
by solid lines. The IRAC and MIPS 24 /irn data points are not connected since the IRAC 
and MIPS data were not taken at the same time within each epoch. 

Figures [7] [8] and El show typical SEDs of evolved stars. The SEDs of O-rich and C-rich 
stars are very similar, both in the shape and in the amount of flux variations observed in 
the IRAC and MIPS 24 /irn bands. The O-rich AGBs tend to have brighter 24 //m fluxes, 
possibly due to the presence of the broad 18 /xrn silicate band. The SEDs of the extreme 
AGB stars tend to peak at longer wavelengths. Extreme AGBs are very bright at 24 /irn 
and some of them show very large variability in all IRAC and MIPS bands. 

Cepheid stars present typical SEDs of stellar photospheres with no evidence for an 
infrared excess (Figure [TUI) . Among the massive stars in our sample, early type stars similarly 
show the characteristic SEDs of stellar photospheres, while RSGs have SEDs similar to the 
ones of O-rich AGBs, with generally small brightness changes in the IRAC and MIPS bands 
(see Figure fill) . The SEDs of YSOs (Figure [T2D are the reddest in our sample, steeply 
increasing from optical to thermal infrared wavelengths. 

Figure [13] shows SEDs of a representative sample of our unclassified sources. Most of 
the unclassified sources appear to be near or on the main sequence in Figure E] panels a and 
b. However, Figure Eb shows that some of these sources appear to be slightly redder (~0.5 
mag in [4.5] — [8.0]) than main sequence stars. It is possible that a fraction of the unclassified 
sources are older YSOs approaching the T Tauri phase, without protostellar envelopes and 
a vestigial disk. Some unclassified sources, however, have SEDs very similar to those of very 
red YSOs, and indeed are found in the same CMD space where variable sources with matches 
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in the YSO catalogs are located. 


5.3. SSTISAGE1C J005212.88-730852.8: a Candidate SMC OH/IR Star 


SSTISAGE1C J005212.88-730852.8 has one of the highest variability indices at 24 /an 
in onr sample: V 24 = —46.25 between Epochs 0 and 1. The negative sign indicates that the 
source has been declining in brightness between November 2004 (S 3 MC) and September 2007 
(SAGE-SMC Epoch 1 ). The variability index in the 24 /mi band between Epochs 0 and 2 
(June 2008) is less negative (V 24 = —30.38), suggesting that the source has been brightening 
in the 9 months between the two SAGE-SMC MIPS 24 /mi epochs. This is confirmed by the 
variability index we measure in all IRAC bands, which is positive between Epochs 1 and 2, 
indicating that the source has continued brightening in the three months between June and 
September 2008 (when the IRAC SAGE-SMC data were collected). 


SSTISAGE1C J005212.88-730852.8 is classified by both Bo yer et ah 120111) and Soszvhsk i et al. 


12 0111) as an O-rich AGB star. The OGLE-III catalog, in particular, classifies the source as 
a Mira with a long secondary period of 1,453 days, favoring its classification as an AGB star 
rather than a RSG (the latter b eing ty pically a Semiregular variable). The O-rich chemistry 
is confirmed bv lvan Loon et al.l (120081) . who note that in the past the star h as be en confused 
with a carbon star. The optical spectral classification by Blanco et al. ljl9801) as an M 6 
giant, however, is in agre ement with the Spitzer and ISO spectra showing a silicate feature 
in emission flCioni et al.l 120031) . The source is peculi ar al so because it is the only stellar 
object in our catalog that has a match in the IWong et ahl (2012) catalog of radio sources: 
the star is detected in the 6 cm band, which is known for the presence of OH masers (the 
source is not detected in the 3 cm band, as expected if the emission is not radio continuum). 
These characteristics suggest that SSTISAGE1C J005212.88-730852.8 may be an OH/IR 
star: if confirmed, this could be the first OH/IR star in the SMC (none are currently known, 
according to van Loon 2012). 


The source is very bright ([3.6] = 8.63 mag and [24] = 4.73 mag), suggesting that it 
may be one of the intermediate mass O-rich AGBs undergoing HBB. It is also very red 
([8.0] — [24] = 2.73 mag), indicative of a large amount of cold circumstellar dust. In fact, 
the SED of this source (see Figure [ 8 ]), is quite similar to the flat, red spectra of extreme 
AGBs. The spectral characteristics suggest that SSTISAGE1C J005212.88-730852.8 may 
be one of the O-rich AGB stars with highest mass loss rate in the SMC, comparable to 
the one typically observed in the C-rich extreme AGB stars. Finding sources of this kind 
is important to explain the overall dust budget in the Magellanic Clouds. As mentioned 
in Section 14.11 the large dust producers in both the LMC and SMC tend to be extreme 
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carbon stars flBover et al.ll2012l; iRiebel et al.ll2012h . making it difficult to explain the average 
composition of the Magellanic Clouds’ ISM (which has similar fraction of O-rich and C-rich 
dust). Stars like SSTISAGE1C J005212.88-730852.8 may help to reduce the observed deficit 
in silicate dust production. 


5.4. SSTISAGE1C J005107.19-734133.3: an RV Tau Variable 


SSTISAGE1C J005107.19-734133.3 is a very red star classified as an extreme AGB by 
Bover et alj (2 0111) and as a Typ e II Cepheid (with a period of 39.52 days) in the OGLE-III 
database (jSoszvriski et ah 2010b). As noted in Section 14721 the star is much redder than all 
other Cepheids, and its SED is dominated by a strong infrared excess indicative of copious 
amounts of circumstellar dust. The absolute variability index of this source at 24 /irn between 
Epochs 1 and 2 is larger than 5, the highest for all sources in the OGLE-III Cepheid class, 
and larger than 7 in all IRAC bands. 


SSTISAGE1C J005 


07.19-734133.3 is in fact among the post-AGB stars spect rally clas- 


sided bv iKamath et al.l (120141) in the SMC. In their analysis, iKamath et al.l (120141) confirms 


that the source shows an RV Tau type variability with alternating shallow and deep minima, 
and an infrared excess consist ent with the presence of a circumbinary disk, as common in 
post-AGB stars. As noted by S oszyriski et ajJ (i2010bi ). RV Tau variables are matched with 
Type II Cepheids in the OGLE-III automatic classification light curve scheme. Visual in¬ 
spection of the optical spectra of this source suggests s-process enrichment, evidence that 
the star has recently gone through the third dredge-up, and thus may have only recently left 
the AGB phase. 


6. Conclusions 

We have identified and characterized a population of candidate infrared variables in 
the SMC using the data from the Spitzer SAGE-SMC and S 3 MC surveys. We found 814 
sources showing variability in two or more adjacent mid-infrared bands. The variable source 
population is dominated by evolved stars (~75%), specifically C-rich AGBs (~60%). Around 
10% of our variable sources are O-rich AGB stars, with one variable star candidate thought 
to be a rare OH/IR source. We End that ~3% and <1% of our variable sources are RSG 
and RGB stars, respectively. Our list of candidate infrared variables include a small, but 
significant population of Cepheid variables (~9%), as well as early type variable stars (~3%), 
and variable YSOs (~5%). A little over 1% of our sources are background galaxies that vary 
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in the infrared. Around 7% of sources in our catalog are currently unclassified, but may be 
further examined through future studies. 
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Table 1. SAGE-SMC and S 3 MC Observing Dates and Survey Areas 


Epoch 

Survey name 

IRAC observation date 

MIPS observation date 

Sky area (deg 2 ) 

0 

S 3 MC 

2005 May 7-9 

2004 November 6-8 

~ 3.7 

1 

SAGE-SMC 

2008 June 12-19 

2007 September 17-25 

~30 

2 

SAGE-SMC 

2008 September 15-23 

2008 June 25-28 

~30 






Table 2. Variability Criteria and Source Counts 


Criteria 

3.6 //in 

4.5 fi m 

5.8 /xm 

8.0 n m 

24 /m i 


Epochs 0 and 1 

# of sources with valid fluxes 

169424 

120190 

27052 

10172 

1425 

# with |v band > 3 

1939 

1708 

570 

409 

189 

% with |v band > 3 

1.14 

1.42 

2.11 

4.03 

13.26 

# of sources with valid fluxes 
in neighboring bands 

115308 a 

116101 b 

23885 c 

11294 d 

1083 e 

# meet variability criterion 

194 

407 

466 

617 

623 

% meet variability criterion 

0.17 

0.35 

1.95 

5.46 

57.53 

Epochs 0 and 2 

# of sources with valid fluxes 

175714 

123027 

27066 

12129 

1551 

# with |v band > 3 

1963 

1742 

571 

421 

189 

% with |v band > 3 

1.12 

1.42 

2.11 

3.47 

12.19 

# of sources with valid fluxes 
in neighboring bands 

112696 a 

113568 b 

23917 c 

9552 d 

976 e 

# meet variability criterion 

196 

411 

483 

645 

646 

% meet variability criterion 

0.17 

0.36 

2.02 

6.75 

66.19 

Epochs 1 and 2 

# of sources with valid fluxes 

812930 

551722 

87886 

36444 

6544 

# with |v band > 3 

3561 

3252 

582 

435 

277 

% with |v band > 3 

0.44 

0.59 

0.66 

1.19 

4.23 

# of sources with valid fluxes 
in neighboring bands 

496470 a 

500512 b 

81995 c 

32302 d 

4092 e 

# meet variability criterion 

446 

446 

390 

374 

227 

% meet variability criterion 

0.09 

0.09 

0.48 

1.16 

5.55 
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“Valid fluxes at 3.6 and 4.5 //m. 
b Valid fluxes at 3.6 and 4.5 /./in or at 4.5 and 5.8 n m. 
c Valid fluxes at 4.5 and 5.8 //in or at 5.8 and 8.0 //m. 
d Valid fluxes at 5.8 and 8.0 //in. 

“Valid fluxes at 8.0 and 24 //rn. 
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Table 3. Properties of the SAGE-SMC Infrared Variables 


Column 

Name 

Description 

1 

varid 

The identification name of the variable star candidates 

2 

R.A. (J2000) 

Right Ascension, J2000 (deg) 

3 

Decl. (J2000) 

Declination, J2000 (deg) 

4-6 

IRACdesignations_e7 

SAGE-SMC IRAC designations for Epoch 0 (7=0), 

Epoch 1 (7=1), and Epoch 2 (7=2) 

7-9 

MIPS designations_e 7 

SAGE-SMC MIPS designations for Epoch 0 (7=0), 

Epoch 1 (7=1), and Epoch 2 (7=2) 

10-12 

IRACdistance_7 

Distance between SAGE-SMC IRAC Epochs 0 and 1 (7=01), 

Epochs 0 and 2 (7=02), Epochs 1 and 2 (7=12) (arcsec) 

13-15 

MIPSdistance_e7 

Distance between SAGE-SMC IRAC and MIPS for Epoch 0 (7=0), 

Epoch 1 (7=1), and Epoch 2 (7=2) 

16 

E01 a 

Variability flag between Epochs 0 and 1 

17 

E02 a 

Variability flag between Epochs 0 and 2 

18 

E12 a 

Variability flag between Epochs 1 and 2 

19-21 

classification) 

Classification of sources (j= 1-3) 

22-26 

V7_ 01 

Variability indices between Epochs 0 and 1 for the 

IRAC 3.6, 4.5, 5.8, 8.0 pm (A;=1-4) and MIPS 24 pm (7=5) bands. 

27-31 

V7_ 02 

Variability indices between Epochs 0 and 2 for the 

IRAC 3.6, 4.5, 5.8, 8.0 pm (A;=1-4) and MIPS 24 pm (7=5) bands. 

32-36 

V7_ 12 

Variability indices between Epochs 1 and 2 for the 

IRAC 3.6, 4.5, 5.8, 8.0 pin (7=1-4) and MIPS 24 pm (7=5) bands. 

37-46 

flux7, dflux7' 

Fluxes and flux uncertainties (mJy) for MCPS UBVI h (7=1-8), 

OGLE-III VI c ' d (7=9-10), 2MASS JHK S (7=11-16) 

47-62 

e0_flux7, e0_dflux7 

Fluxes and flux uncertainties (mJy) for Epoch 0: 

Spitzer IRAC 3.6, 4.5, 5.8, 8.0 pm (7=1-4) and MIPS 24 pm (7=5) bands. 

63-72 

el_flux7, el_dflux7 

Fluxes and flux uncertainties (mJy) for Epoch 1: 

Spitzer IRAC 3.6, 4.5, 5.8, 8.0 pm (7=1-4) and MIPS 24 pm (7=5) bands. 

73-82 

e2_flux7, e2_dfiux7 

Fluxes and flux uncertainties (mJy) for Epoch 2: 

Spitzer IRAC 3.6, 4.5, 5.8, 8.0 pm (7=1-4) and MIPS 24 pm (7=5) bands. 

83-96 

mag k, dmag7 

Magnitudes and magnitude uncertainties (mag) 
for MCPS UBVI (7=1-8) and 2MASS JHK S (7=9-14) 

97-106 

e0_mag7, e0_dmag7 

Magnitudes and magnitude uncertainties (mag) for Epoch 0: 

Spitzer IRAC 3.6, 4.5, 5.8, 8.0 pm (7=1-4) and MIPS 24 pm (7=5) bands. 
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Table 3—Continued 


Column 

Name 

Description 

107-116 

el_magZ, el_dmagZ 

Magnitudes and magnitude uncertainties (mag) for Epoch 1: 

Spitzer IRAC 3.6, 4.5, 5.8, 8.0 //in (k= 1-4) and MIPS 24 //m (/=5) bands 

117-126 

e2_mag k, e2_dmagZ 

Magnitudes and magnitude uncertainties (mag) for Epoch 2: 

Spitzer IRAC 3.6, 4.5, 5.8, 8.0 //m (k= 1-4) and MIPS 24 //in (/=5) bands 


a Flags indicating if a given source has been identified as a candidate infrared variable based on 
the data from two epochs (E01 - Epochs 0 and 1; E02 - Epochs 0 and 2; E12 - Epochs 1 and 2): 1 
- a source is a variable, i.e., its absolute variability index is larger than 3 in at least 2 consecutive 
bands; 0 - a source does not vary or there is not enough photometric data available. 

b The MCPS catalog provides magnitudes only. The U-, />-, V-, and /-band fluxes and flux 
uncertainties were derived using the zero magnitude flux of 1790, 4063, 3636, and 2416 Jy, respec¬ 
tively. 

c No photometric uncertainties are available for the OGLE-III data. 

d The OGLE-III catalog provides magnitudes only. The V- and /-band fluxes were derived using 
the zero magnitude flux of 3636 Jy and 2416 Jy, respectively. 
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Table 4. Catalog Matching Results 


Literature 


Matched contents # of sources 


Bolatto et al. (2007) 

Bonanos et al. (2010 1 

Bover et al. (2011 1 

Carlson et al. (2011 1 

G. Jacoby (private communication) 

Kozlowski et al. (2011 1 

Kozlowski et al. (2013 1 

Oliveira et al. (2011 1 

Oliveira et al. (2013 1 

Seale et al. (20141 

Sewilo et al. (20131 

Sheets et al. (20131 

Simon et al. (20071 

Soszvnski et al. (2010a) 

Soszvnski et al. (2010b 1 
Soszvnski et al. (2011 1 
Sturm et al. (20131 
van Loon et al. (2008 1 
van Loon et al. (2010 1 
Wong et al. (20121 


YSOs 

13 

Massive stars 8. 

43 

Evolved stars b,f 

600 

YSOs 

3 

PNe 

1 

AGN C 

8 

QSOs c 

1 

YSOs 

1 

YSOs 

4 

YSOs 

10 

YSOs 

44 

dustyOB d 

1 

YSOs 

3 

Classic Cepheids 6 

67 

Type II Cepheids 6 

3 

Evolved stars f 

539 

X-ray sources 

6 

YSOs 

1 

YSOs 

1 

Radio 

1 


includes early type and RSG stars. 

b Includes cAGB, xAGB, oAGB, RGB and RSG stars. 

c Reported as a background galaxy. 

d Reported as an early type star. 

e Reported as a Cepheid. 

includes cAGB, xAGB, oAGB and RGB stars. 
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Fig. 1.— The three-color composite image of the SMC combining the SAGE-SMC Epoch 
1 MIPS 24 pm (red), IRAC 4.5 prn (green), and IRAC 3.6 pm (bine) images. The white 
contour outlines the S 3 MC survey area. The main components of the SMC (Bar, Wing, and 
Tail) are indicated. 
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Fig. 2.— Distribution of the absolute values of variability indices (|Vj) in each IRAC and 
MIPS 24 /ini band, calculated for the interval between Epochs 0 and 1 (top), 0 and 2 
(middle), and 1 and 2 (bottom). Measurements between Epochs 0 and 2 yielded a similar 
spread of |V| as between Epochs 0 and 1; both sample variability over a period of a few 
years. Measurements between Epochs 1 and 2 instead sample a change over a period of 
just a few months. Sources with |V| > 3 in at least two consecutive bands are classified as 
variable sources. 
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Fig. 3.— Distribution of the absolute values of variability indices (|1/|) between Epochs 1 
and 2 for the 1RAC 3.6 /irn band for several classes of sources as indicated in the legends. 
The largest variability indices are found for Cephcids and extreme AGB stars. The vertical 
line corresponds to |E|=3. 
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Fig. 4.— Median value of the absolute variability index (|G|) for each band, each pair of three 
epochs, and each source class. The standard deviations of the variability index distributions 
are represented by vertical bars. The variability indices are larger for the longer inter-epoch 
interval, sampling time-scales of ~3 years, with the exception of Cepheids that have similar 
variability in all intervals. O-rich AGBs and YSOs have the largest spread in |V| at 24 /jrn, 
while C-rich AGBs have the largest spread in all bands. 
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Fig. 5.— The spatial distribution of the variable star candidates in the SMC. The sources 
are overlaid on the SAGE-SMC 8.0 pm image. Sources with multiple matches have all 
classifications plotted. 







39 



5 6 


-2 


J-K* 


4 6 

J-I3.6] 


10 



12 3 

[4.5M8.0] 


4 5 


-1 


12 3 4 

[8.0H24] 


Fig. 6.— The K s vs. J—K s (a), [3.6] vs. J-[3.6] (b), [8.0] vs. [4.5]-[8.0] (c), and [8.0] vs. [8.0]- 
[24] (d) CMDs showing the distribution of the variable star candidates with respect to the 
entire population of sources from the SAGE-SMC IRAC Epoch 1 catalog (shown in grey as 


Sewilo et al. 

2013; 

Chen et al. 

2014) 


evolved stars (e.g.. iBover et al.l 20111) and YSOs (e.g. 
respectively. The symbols represent different classes of objects as indicated in the legends 
Sources with multiple matches have all classifications plotted. 
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Fig. 7.— SEDs showing a representative sample of variable sources classified as C-rich AGB 
stars. The SAGE-SMC IRAC and MIPS 24 /jrn Epoch 0, Epoch 1, and Epoch 2 fluxes 
are plotted as red squares, blue circles, and pink triangles, respectively. The ancillary data 
from MCPS, OGLE-III, and 2MASS/6X2MASS surveys (see Section 2.1) are shown as black 
circles. The data obtained with the same instrument are connected by solid lines. 
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Fig. 8.— SEDs showing a representative sample of variable sources classified as O-rich AGB 
stars. The symbols and lines are the same as in Figure [0 
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Fig. 9.— SEDs showing a representative sample of variable sources classified as extreme 
AGB stars. Most of these sources have a C-rich chemistry. The symbols and lines are the 
same as in Figured 
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Fig. 10.— SEDs showing a representative sample of sources classified as Cepheid variable 
stars. Note that the very red source SSTISAGE1C J005107.19-734133.3 is not a Cepheid, 
but rather a post-AGB star with RV Tau variability (see Section 15.41) . The symbols and 
lines are the same as in Figure 0 
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Fig. 11.— SEDs showing a representative sample of variable sources classified as massive 
stars. SEDs of early type stars and RSGs are shown in the left and right panel, respectively. 
The symbols and lines are the same as in Figure [71 
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Fig. 12.— SEDs showing a representative sample of variable sources classified as YSOs. The 
symbols and lines are the same as in Figure [71 
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Fig. 13.— SEDs showing a representative sample of unclassified variable sources. The 
symbols and lines are the same as in Figure [71 


































